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ABSTRACT 

The power spectrum of cosmic infrared background (CIB) anisotropies is sensitive to 
the connection between star formation and dark matter halos over the entire cosmic 
star formation history. Here we develop a model that associates star-forming galaxies 
with dark matter halos and their subhalos. The model is based on a parameterized 
relation between the dust-processed infrared luminosity and (sub)halo mass. By ad- 
justing 3 free parameters, we attempt to simultaneously fit the 4 frequency bands of 
the Planck measurement of the CIB anisotropy power spectrum. To fit the data, we 
find that the star-formation efficiency must peak on a halo mass scale of « 5 x 10^^ 
and the infrared luminosity per unit mass must increase rapidly with redshift. By com- 
paring our predictions with a well-calibrated phenomenological model for shot noise, 
and with a direct observation of source counts, we show that the mean duty cycle of 
the underlying infrared sources must be near unity, indicating that the CIB is dom- 
inated by long-lived quiescent star formation, rather than intermittent short "star 
bursts" . Despite the improved flexibility of our model, the best simultaneous fit to all 
four Planck channels remains relatively poor. We discuss possible further extensions 
to alleviate the remaining tension with the data. Our model presents a theoretical 
framework for a future joint analysis of both background anisotropy and source count 
measurements . 

Key words: submillimetre: diffuse background - submillimetre: galaxies - galaxies: 
star formation - galaxies: halos 



1 INTRODUCTION 



About half of the optical/UV emission from stars is absorbed 
by dust in galaxies and re-emitted a t infrared wavelengths 
l|Dwek et al.ll 19981 : iFixsen et al.lll998t ). Measurements of the 
CIB therefore offer an important channel to study the star- 
formation history of the Universe. In particular, CIB mea- 
surements are sensitive to the star formation out to high 
redshifts, due to the correlation between the infr ared lumi- 
nosit y of a galaxy and its star formation rate l|Kennicuttl 
Il998f ). Owing to the relatively poor angular resolution of 
far-infrared telescopes, it is challenging to resolve and study 
individual infrared sources. Surveys to date are only able 
to resolve the brightest sources, which are responsib le for 
less than 15% of the total CIB (|01iver et all l2010bl ). Al- 
though the Atacama Large Millimeter/Submillimeter Array 



(ALMAfl will soon be able to resolve most of the back- 
ground, it will only do so over a small region of the sky. 
Large-scale fluctuations of the unresolved background will 
contain complementary information that will remain use- 
ful even after the advent of ALMA (|Haiman fc Knoxll2000l : 
iKnox et al]|200ll ). 

Althou gh fluctuations in the CIB were detected over a 
decade ago (|Lagache fc Pugedl20odl : iMatsuhara et aLlbOOOD . 
these were limited to probing only the Poisson contribution 
to the power spectrum. Robust measurements of the clus- 
tered component in the CIB anisotropy, probing the spatial 
correlations of the underlying sources, have become possible 
relatively recently, and are now rapidly improving. The first 
detection of clustered power came from 160 ^, T[^ Spitzer data 
(IGrossan fc Smootll2007l : iLagache et alll2007l ) . This was fol- 
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lowed recently by the Balloon-borne Large Aperture Submil- 
limeter Telescope (BLASTfl which measured the clustering 
power near the energetic peak of the CIB at 6 00, 860, and 
1200 GHz (500, 350, an d 250 /im, respectively: IViero et all 
l2009l ). lHall et all l|2010l) reported the first detection of clus- 
tered CIB power at millimeter wavelengths, in South Pole 
Telesc ope (SPTjl data, soon followed by iDunklev et alj 
('2010') i n Atacama Co smology Telescope (ACT^f data, and 
[Shirokoff et al.l l|201ll ') with additional SPT data. Most re- 
cently, measurements of the CIB fluctuation power, extend- 
ing over a broad ran ge in frequencies and angular scales, 
have be en reported in lAde et al.l (|201ll ) and I Amblard et "all 
l|2011a|) . inferred from Planck and Herschel data, respec- 
tivelylfl In addition, a joint analysi s of multiple measure - 
ments has also been performed by lAddison et al.l (|201lD : 
they find that the power spectra over the range of angu- 
lar scales and frequencies considered could be well fitted by 
a simple power law. 

To make full use of this information, and to draw conclu- 
sions about the cosmic star formation history, requires the 
development of a modeling framework. With too complex a 
model, degeneracies between model parameters threaten to 
make meaningful conclusions impossible. On the other hand, 
if the model is too simple, it may be unable to accommo- 
date all relevant information, such as observations of bright 
individual sources - or even worse: we may risk having con- 
clusions that are artificially driven by inadequate modeling 
assumptions. 

The Planck and Herschel measurements, the most infor- 
mative to date, have so far been interpreted vi a halo models 
of large-scale structure jCoorav fc Shethll2002l ) . These mod- 
els take advantage of the relatively small uncertainty with 
which, given a model, statistical properties of the spatial 
distribution of dark matter can be calculated. They do so 
by dividing the modeling into two stages. First, the spatial 
distribution of dark matter is described (in particular, by 
assuming that all dark matter is confined to lie with spher- 
ically symmetric, collapsed halos). Second, the observable 
of interest is related to the distribution of dark matter. A 
key element for the latter is the halo occupation distribution 
(HOD) function of galaxies, which is a statistical description 
of how galaxies occupy halos, depending on halo mass and 
redshift. 

To avoid degeneracies that might arise in th e presence 
of too ma ny free paramet e rs, the halo model in lAde et ahl 
l|201ll ) and lAmblard et al.l (|2011al '). to which we refer as the 
"standard model" in the following, has been relatively sim- 
ple - ignoring many details of how infrared-luminous galax- 
ies trace dark matter halos. For example, the HODs were 
assumed to be the same for all spectral types and to not 
vary with redshift. Even quite coarse details, such as the 
dependence of the luminosity of a galaxy on halo mass have 
been neglected. 

These models have the virtue of simplicity, and without 
the data demanding more complicated models, it would per- 
haps not be worth abandoning such simplicity. The primary 
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motivation of our study, however, is that the existing data 
are already difficult to understand with these simple models. 
For example, in the models used to interpret the Planck CIB 
measurements, there is no single choice of parameters that 
generates a goo d fit, simultaneo usly, to the power spectra at 
all frequencies l|Ade et al.ll201ll ). Also, because these models 
do not include a dependence of the IR luminosity on halo 
mass, neither the number counts of individual sources (de- 
tectable at the bright end), nor the shot-noise power levels 
can be calculated reliably. 

In this paper, we extend the modeling framework used 
to date, by allowing galaxies in halos with different masses 
to have different luminosities. We do so by assuming that 
satellite galaxies occupy dark matter subhalos, and use a 
subhalo mass function derived from cosmological N-body 
simulations. The luminosity of each satellite galaxy is then 
parametrically related to the mass of the subhalo. Ours is 
not the first CIB mo del to allow for IR galaxi es to have 
different luminosities. lAmblard fc Cooravl (|2007l ) discussed 
such a model, based on t he co n ditional luminosity function 
(Coorav & Milosavlicvi 3 I2OO5I : lYang. Mo fc van den BoschI 
2003. : ,Yang et a l. 2005). What is new in the present work 
is the explicit association of IR galaxies with subhalos, and 
the use of such a model for interpreting the Planck data. 
As we will show below, the comparison of this model to the 
data indeed leads to novel conclusions. 

The modeling extension necessarily comes at the cost of 
a few additional parameters, but has several benefits. First, 
and perhaps most importantly, it allows us to estimate a 
lower limit to the bright end of the number counts. Sec- 
ond, by introducing a duty cycle, and comparing to source 
count observations, we can put limits on the duty cycle. In 
the context of our model, the data favor duty cycles near 
unity. Physically, this implies that the bulk of the CIB is 
produced by quiescent star formation, rather than by the 
intense, sho rt starbursts expec t ed to be triggered in ma - 
jor mergers (|Sanders et al.|[l98^ : lBarnes fc Hernquistlll99ll '). 
Third, there exists a tension between the value of the power- 
law index a, describing the dependence of the number of 
satellite galaxies on halo mass, inferred either from cosmo- 
logical simulations, from optical galaxy data, or from the 
CIB anisotropy measurements (see Our model naturally 
resolves this tension. 

The remainder of this paper is organized as follows. In 
321 we review recent measurements of the CIB anisotropy 
power spectrum by BLAST, SPT, and the space missions 
Herschel and Planck. In !j3l we first review the basic for- 
mulae of the existing halo-based CIB models, and then de- 
scribe our extensions that include a more realistic L — M 
relation. In ^ we discuss the ingredients of the extended 
model in more detail, and illustrate how each ingredient af- 
fects the resulting CIB anisotropy power spectra. In fjH we 
apply our model to the recent measurement by Planck, and 
derive constrains on both the L — M relation and on the 
shape of the mean spectral energy distribution (SED) of the 
sources. With enhanced flexibility, we attempt to simulta- 
neously fit the data at all measured frequencies. In IJSl we 
examine the shot noise levels and source counts computed 
with the new model, compare them with an observationally- 
calibrated model and a direct measurement, and place a 
constraint on the duty cycle of the underlying infrared 
sources. In iJTl we further discuss the details of the fit, as 
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well as several caveats and potential future improvements 
to the model. Finally, in ^ we summarize our main conclu- 
sions. Throughout this paper, we adopt the standard ACDM 
model as our fiducial background cosmology, with parame- 
ter values consistent with the seventh year WMAP results 
l|Komatsu et al.ll201ll ). i.e., {nm,QA,0,t,h,a8,n.s} = {0.274, 
0.726, 0.045, 0.705, 0.810, 0.96}. 



2 MEASUREMENTS OF CIB ANISOTROPIES 

Substantial progress has been made in measuring the power 
spectrum of the CIB anisotropies in the past few years, 
thanks to the joint efforts from BLAST, SPT, Herschel and 
Planck. 

In 2009, BLAST measured the power spectra in the 
GOODS-South field at 250, 350, and 500 fim (i.e., 1200, 857, 
600 GHz in frequencie s, respectively) in the multipole range 
of 940 ^ Z < 10,800 l|Viero et al.H2009l '). After subtracting 
the Galactic cirrus and Poisson noise, they found the vari- 
ance of the CIB over the scales of 5' - 25' is consistent with a 
constant amplitude of 15% ± 1.7% relative to the CIB mean 
at all observing frequencies. The data could alternatively be 
well fit by a linear halo bias model, with bias parameters, 
6 = 3.8 ± 0.6, 3.9 ± 0.6 and 4.4 ± 0.7 at 1200, 857, 600 GHz, 
respectively. Interestingly, the data could not be fit by the 
standard halo model based formula, with data points lying 
below the model curves at small scales and above the model 
curves at large scales. Later comparison with the Planck 
measurement shows a general consistency between the two 
measurements. However, the data points from the BLAST 
measurement are systematically higher at large scales. This 
seems to suggest that the BLAST measurement have been 
contaminated by residual Galactic cirrus at the large scales 
l|Ade et al.ll2oTlh . If true, the failure of the halo model fit, 
likely due to the contamination of the data, does not neces- 
sari ly mean the model is wrong. 

Hall et al.1 (|2010l ) reported SPT measurements of the 



auto- and cross-correlation at 150 GHz and 220 GHz in the 
multipole range of 2000 ^ I ^ 10000. They also determined 
that the spectral indices of the Poisson and clustered compo- 
nents between the two observing frequencies are 3.86 ± 0.23 
and 3.8 ± 1.3, respectively. This implies a steep slope of the 
SED of the dust emission (i.e. a gray-body power-law index 
of /3 ~ 2; see definition in eq. [24] below). Although the SPT 
measures the small-scale clustering, a model based on linear 
halo bias was still able to provide an acceptable fit to the 
data. This does not necessarily mean the clustering is lin- 
ear at such small scales, but rather shows that the shape of 
non-linear power spectrum over the scales probed by SPT 
can not be distinguished from that of the linear power spec- 
trum. A clear detection of the excess power from non-linear 
clustering requires the measurement to cover the transition 
at multiples of I ~ 1000. This is exactly the angular coverage 

provided by Herschel and Planck. 

Soon afterwards, iDunklev et al.l (|2010l ) published their 
measurement of the power spectra at 148 GHz and 218 GHz, 
based on the ACT data collected during their 2008 season 
over 296 deg^. The clustered component of the IR power 
spectrum is detected at 5a, assuming an analytic model for 
its shape. The spectral index of unresolved IR emission be- 



tween the two frequencies is found to be 3.69 ± 0.14, consis- 
tent with the SPT measurement. 

With the Spectral and Photometric Imaging Receiver 
(SPIRE) on board, the Herschel Space Observatory, during 
its Science Demonstration Phase, measured the power spec- 
tra of CIB anisotropies in the Lockman Hole and GOODS 
South field at 250, 350, and 500 fj ,m over the scales 1' - 
40' (540 ^l ^21600 in multipoles: lAmblard et ai] l2011bD . 
The coverage of a wide range of angular scales, combined 
with clean data, yielded a clear detection of the non-linear 
clustering signal. The standard halo model-based prescrip- 
tion provides satisfactory fits to the data, but with different 
sets of best-fit parameters at different frequencies. The min- 
imum halo mass for star-formation, M,nin , is constrained to 
be ~ 3 X IO'^^Mq, which the authors interpret as the most 
efficient mass-scale for star formation. This mass-scale is 
lower than previous predictions by semi-a nalytical models 
for galaxy formation ([Gonzalez et al.l|2011^ . 

Most recently, the Planck team published their mea- 
surement of the angular power spectra of CIB anisotropies, 
using maps of six regions of low Galactic dust emission 
(Ade et al. 20lJ, hereafter All). The power spectra were 
determined in 4 frequency bands (217, 353, 545, and 857 
GHz) over multipoles between / = 200 and I = 2000. Planck 
found CIB anisotropies at the level AI/I — 15% com- 
pared to the mean CIB, and also that the fluctuating com- 
ponent and the mean CIB have similar frequency spectra, 
consistent with the BLAST results. Unlike other missions, 
Planck was unable to independently measure the amplitude 
of the shot noise, due to its relatively poor angular reso- 
lution (~ 5 arcmins). Instead, it adopted values calc ulated 
using a parametric model by (|Bethermin et al]|201ll . here- 
after Bll). Again, different frequency bands yield different 
best-fit model parameters. 



3 HALO MODEL FOR CIB ANISOTROPIES 
AND ITS EXTENSION 

3.1 Review of the existing models 

The angular power spectrum of CIB anisotropies is defined 
through. 



(1) 



where f denotes the observing frequency and Ii, is the mea- 
sured specific intensity at that frequency. In a flat universe 
as assumed throughout this paper, the specific intensity is 
related to the comoving specific emission coefficient j via 



dz •'^ ' ' 



(2) 



d4a,>,.) (l + ^ 
dz \ 3{v,z) 



where x[z) is the comoving distance to redshift z, and a = 
1/(1-1-2:) is the scale factor. Combining eq uations ffl and © 
and employing the Limber approximation (|Limbeil 19541 ) . we 
obtain 



C, 



'^'^a^3{v,z)j{v' ,z)Pj^^^,{k = l/x,z), (3) 



4 Shang et al. 



where is the 3-D power spectrum of the emission co- 

efScient, and is defined as follows, 

< Sj{k,u)Sj{k' ,1^') >= {2T,f]{j.)]{i^')P,^,,,{k)5\k~k'). (4) 

The existing models equate Pj with the galaxy power spec- 
trum Pgai, assuming the CIB is sourced by galaxies, and 
that the spatial variations in the emission coefficient trace 
the galaxy number density, 



(5) 



On large scales, in the linear regime, the galaxy power spec- 
trum follows the linear matter power spectrum, Pgai = 
b^Piin, while on small, non-linear scales, Pgai is typically 
computed from the halo model. In the framework of the 
halo model, the galaxy power spectrum is a sum of three 
terms 



Pgalik, z) = Plh{k, z) + P2h{k, Z) + Pshot{k, z), 



(6) 



where Pih and P2h account for contributions from galaxy 
pairs in the same halo and in different halos, respectively, 
and Pshot is the shot noise. The analyt ical expressions for 
Pih and P2h are jCoorav fc Shethll2002h . 



PiH{k,z) = J dM^{M, z) X 
2N,,„{M)Nsat{M)u{k,z\M) + Nlt{M)u\k,z\M) 



(7) 



'■gal 



and 



P2h{k,z) = Pm{M,z) X 
dM [M, z) ^1^^^ ' " 

dM ' flgal 



(8) 



where M is the halo mass, dN/dM is the halo mass function, 
u{k, z\M) is the Fourier transform of the halo density profile, 
and Ncen and Nsat, specified by the HOD, are the number 
of central and satellite galaxies inside a hal o, with Nggi = 
Ncen+ Nsat- Motivated by simulations (e.g.. iKravtsov et al.l 
I2OO4I '). Ncen is typically modeled as a simple step function. 



M < Mmin 



while Nsat is parameterized by a power law, 
M 



Nsat 



Msat 



(9) 



(10) 



Here, Msat is the arbitrary "pivot" halo mass that hosts 
one satellite galaxy. Note that in the above equations, Nsat, 
Ncen, and Ngai could be thought of as the average number 
of galaxies inside halos with a fixed mass M, with negligible 
scatter. The parameterizations above have been extended to 
include non-zero stochasticity in the relations hip between 
halo mass and the number of galaxies (e.g., IZheng et al.l 
l2005h . 



3.2 The improved model 

The model above assumes that emissivity density traces 
galaxy number density (equation [S| . This assumption im- 
plies that all galaxies contribute equally to the emissivity 
density, irrespective of the masses of their host halos. In 



other words, it assumes that all galaxies have the same lumi- 
nosity. Realistically, however, both the luminosity and the 
clustering strength are closely related to the mass of the 
host halo. In particular, galaxies in massive halos are likely 
both luminous and highly clustered. In general, introducing 
a monotonic L — M relation should give rise to stronger clus- 
tering than the existing models, in which the more strongly 
clustered sources (massive galaxies) have been assigned ar- 
tificially low luminosities. This improvement is indeed the 
main new feature of our model. 

In the following, we abandon the assumption of a mass- 
independent luminosity, and use the standard halo model to 
directly compute the power spectrum of the emission coeffi- 
cient, Pj. Note that the emission coefficient is related to the 
underlying galaxy population as follows. 



(11) 



where L denotes the infrared luminosity and dn/dL is the 
infrared galaxy luminosity function. Neglecting any scatter, 
the galaxy luminosity is a function of the mass of the host 
dark matter halo (for central galaxies) , or subhalo (for satel- 
lite galaxies), and equation (lllf) can be re- written as: 

>(2) = J dM^{z)^[NcenLcen,(l + .)AM,z) (12) 

+ j dm^{M,z)Lsat,[i+z)i.{'m) , 

where m is the subhalo mass, and dn/dm is the subhalo 
mass function. Studies show that the luminosity of a satel- 
lite galaxy best correlates with the mass or circular ve- 
locity of the host subhalo at the time it is accreted into 
the main halo, i. e. before it looses significant mass due 
to tidal s tripping jNagai fc Kravtsov 20051: Vale fc OstrikeJ 
|2006; Co nrov. Wechsler fc Kravtsovll2006l : IWang et al.ll2006l : 
Wet zel fc Whitell201 0'). We therefore use this "unstripped" 
mass to infer the infrared luminosity. 

The computation of the power spectrum of the emission 
coefficient closely follows that of the galaxy power spectrum 
in the previous section. The only difference is that the galaxy 
numbers {Ncen and Nsat) should be replaced by expressions 
accounting for contributions from galaxies to the overall 
emissivity. Note that apart from the additional weighting 
by luminosity, equation (|12|) is essentially identical to the 
expression for the number density: 



n{z) 



dM^{z) [Ncen (M, Z) + Nsat (M, z)]. (13) 



To compute Pj , we therefore only need to replace Ncen by 



fnM,z) = Nc, 



L(i+,),{M,z) 



47r 



(14) 



and Nsat by 
fT\M,z)= dm 



dUsubhalo , .,,-.L(i + z)v(m,z) 

(m,z\M)^ 7- (15) 



dm 



4:TT 
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in equations (O and ((Sj. The final results are 



Plh,vv' {k, z) 



dM 



dN 
dM 



(16) 



P2h,vv' {k, z) 



where 



x[fr{M,z)rj\M,z)u{k,M,z) 
+f:r{M,z)r\M,z)u{k,M,z) 
+fT\M, z)rj\M, z)u'{k, M, z)], 

J^D,{k,z)D,,{k,z)Pu„{k,z), (17) 



.dN . 



D4k,z) = / dM—b{M,z)u(k,M,z)) (18) 

x[/r"(M,z)+/r*(M,z)]. 

It is straightforward to show that equations (1161) and (jlTp 
reduce to equations ((TJ and (O for a flat L — Ad relation 
(galaxies of different masses have the same luminosity), as 
assumed in previous models. 

The main ingredients of this model, such as the 
halo mass function, halo bias, and halo density profiles, 
have been carefully calibrated using numerical simula- 
tions. Throughout this study, we define halos as over- 
dense regions with a mean density equal to 200 times the 
mean density of the universe. W e assume an NFW profile 
l|Navarro. Frenk fc White! 1 19971) for the halo density, and 
adopt the fitting function of Tinker et al.l (|2008l ) for the halo 
mass function and its associated pre scription for the halo 
bias l|Tinker. Wechsler fc Zhenj2010l '). For the subhalo mass 
functi on, we use the fitting function of [Tinker fc Wetzel 
l|201(]| . equation (12) in their paper). 

We note that a luminosity-weight ing scheme s imilar 
to the abov e have been explored by IShethI (|2005l ) and 
ISkibba et al.l (2006), and applied to study resolved sources 
and the environment-dependence of galaxy properties. To 
the best of our knowledge, this is the first time that such a 
weighting scheme has been used for an analysis of unresolved 
brightness fluctuations. 

Finally, in addition to the clustering caused by corre- 
lated large-scale structures, measurements of the CIB power 
spectrum include shot noise from random fluctuations in the 
discrete number of galaxies. In principle, this shot noise can 
be computed in the above model, through the equation 



^shot 



(19) 



where A'^ includes both halos and subhalos, and S is the 
source flux. 



aL,i 



(l + z)i. 



(20) 



and Scut is the flux above which individual sources are de- 
tected and removed in a given experiment. Unfortunately, 
this computation will likely remain inaccurate, given the 
simplicity of our model. In particular, our model neglects 
any scatter in the L — M relation, for simplicity. This scatter 
is particularly relevant to the calculation of the shot noise: 
as is clear from equation (|19[) . shot noise involves the inte- 
gration of L^, which increases with scatter for a fixed mean 
L — M relation. By neglecting scatter, our model therefore 
underestimates the shot noise, and our result should be un- 
derstood as a lower limit. For this reason, we adopt the 



values computed using the parametric model of Bll (their 
model successfully fits bright source count measurements, 
and is also adopted by All as the shot noise) when we fit the 
measured angular power spectra. The values can be found 
in in Table ^ below (and also in Table 3 of All). How- 
ever, later in § [S] equation (|19p and (|20p are employed to 
constrain the duty cycle of the underlying sources. 



4 MODELING THE L M RELATION AND 
ITS EFFECT ON CIB ANISOTROPIES 

4.1 Parameterizing the L ~ M relation 

Before applying our model to the Planck data, we here spec- 
ify and discuss the model ingredients in detail. In the halo 
model, the galaxy power spectrum is fully determined by the 
HOD, namely the functions Ncm(M,z) and Nsat{M, z). In 
our model for the background fluctuations, the power spec- 
trum depends, additionally, on the function Li^i^^^^{m, z). 
The latter depends on three variables: the redshift z, the 
mass of the host (sub)halo, and the observing frequency f. 
Limited by the current data quality, we adopt a relatively 
simple form for L(i+2)^(m,z) in this study. First, we do 
not make a distinction between halos and subhalos of the 
same mass, i.e., we assume that the L — M and L — m le- 
lations are identical. Second, we assume, for simplicity, that 
all galaxies have the same SED, independent of their masses 
and redshifts, and that the L — M relation does not evolve 
with redshift except for an overall normalization. The de- 
pendence of L on the three variables is then separable, and 
L(i+z)v{Tn, z) can be written as the product 



LoHz)^im)e[{l + z)u], 



(21) 



where Lo is a normalization factor. In the following, we dis- 
cuss each of the three components in more detail. 

(1) Redshift evolution $(z) 

For a given (sub)halo mass, the luminosity and the star 
formation rate (SFR) are expected to increase with redshift, 
because of the higher gas accretion rates, higher gas frac- 
tions, and more compact geometries at early times. Addi- 
tionally, mergers, which can trigger starbursts, are more fre- 
quent at high redshift. In this study, we parameterize ^{z) 
as a power law. 



<i>{z) = (14-^) 



(22) 



This is partly motivated by the study of the specific star 
formation rate (sSFR; SFR per unit stellar mass). If we as- 
sume the stellar mass to halo mass ratio evolves only mildly 
with redshift (supp orted by semi-analytical studies such as 
iNeistein et al.|[201ll ). the sSFR should have a redshift evolu- 
tion similar to that of Lia/Af due to a correla tion between 
SFR and infrared luminosity (|Kennicuttlll998t ) . 

The evidence for a smooth power-law evolution 
of the sSFR with redshift is somewhat mixed. Semi- 
analytical models of galaxy formation indeed show that 
the redshift evolution of sSFR, assumed to be primar- 
ily driven by the fresh gas supply, cou ld be well fit by 
a po wer-law with a slope o f ~ 2 .5 ([Neistein fc Dekell 



l2008l : iDekel. Sari fc Ceverinol l2009l : lOliver et al 

Observations, however, suggest a more complicated 



2010al ). 
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shape, with a steep evolution bel ow z <2, and a 
plat eau beyond this redshift (see Bouche et al.l I2OI0I 
and IWeinmann. Neistein fc Dekell I 20 111 for a compilation 
of measurements). Such a transitio n is hard to ex- 
plain from a theoretical perspecti ve (| Bouche et al.l I2OI0I : 
IWeinmann. Neistein fc Dekell l201lh . Due to the different 
shapes, theoretical models tend to predict a sSFR lower by 
a factor of a few at z ~ 2 and much higher at z > 4 than 
inferred from observations, even though these models could 
more or less reproduce the redshift evoluti on of the over- 
aU SFR (|Weinmann. Neistein fc Dekell [ioTH '). On the other 
hand, there are large uncertainties in current observations, 
especially a t high reds hifts. For example, measurem ents by 
lYabe et al.l Hooi) and lSchaerer fc de Barro^ (|2010l ') suggest 
the sSFR continue to rise beyond z > 2. Given these uncer- 
tainties, we consider two scenarios: Sz has a single value at 
all redshifts (case 0), and is fixed to at 2 > 2 (case 1). 

(2) L-M relation E(m) 

A robust conclusion from comparisons between ob- 
served galax;y luminosity functions and halo mass func- 
tions is that star formation is only effective in a cer- 
tain range of halo masses. In particular, star-formation 
must be suppressed at both low and high halo masses - 
for example, by feedback processes such as photoioniza- 
tion heating, supcrnovac, active galactic nuclei (AGN) and 
virial shocks ( Birnboini fc Dekel 200p: [Keres et all I2OO5I: 



Dekel fc BirnboimI bood : iBower et al.ri2006l : ICroton et all 
20061 ). The quantity of particular interest is the halo mass 
scale at which the star formation is most efficient, corre- 
sponding to a peak in the luminosity-to-mass (L/m) ratio. 
Assuming, for simplicity, that L/m as a function of mass m 
is log-normal, we model the L — m relation as follows, 



E(m) 



27ra? , 

L/m 



: exp 



(logio m - logio m^fff 



2a? , 

L/n 



(23) 



Here rrieff describes the peak of the specific IR emissivity 
per unit mass, and cr^ describes the range of halo masses 
for producing IR luminosity. At the low-mass end, in ad- 
dition to the exponential drop in equation (|23|) . we impose 
a minimum mass of IO^^Mq, i.e., L = for m < 10^^M(7) . 
This cutoff is motivated by the work of lBouche et al.l (|201Cf ) . 
Later, we will see the best-fit rrii-f f is more than 2(7^ away 
from the cutoff, so effects of applying such a cutoff are mild. 

(3) SED shape e{iy) 

W e adopt the same form for the SED as iHall et all 
l|2010l ). The SED is a graybody at low frequencies, and is 
a power law at high frequencies. 



(24) 



where is the Planck function and Td is an effective dust 
temperature. These two functions are connected smoothly 
at the frequency that satisfies 



= -7. 



(25) 



We fix 7 = 2 as in iHall et all (|2010l '). but, unless stated 
otherwise, we allow /? and Td to be free parameters. 



4.2 Effects on CIB anisotropics 

Next, we illustrate how the CIB anisotropy power spectrum 
depends on each of the model parameters. For this purpose, 
we adopt the following set of fiducial model parameter val- 
ues: = 3, me// = 10'2 '5 M0, ai/^ = 0.5, Ta = 34 K, 
P = 2. The choice of = 3 is based on low redshift measure- 



ments of the sSFR (iNoeske et al.ll2007l: iDunne et all [20091 : 
lOliver et alll2010al : iRodighiero et al.l l2010h : while m./ / = 
IqI^-S jg around the peak of the stellar-mass to halo 
mass ratio found by the abundance matching technique 
(e.g., IGuo et al.l [2010l ). The val ue Tg = 34 K is consis- 
tent with measurements such as iDunne et al ] (|200d . Td = 
35.6 ± 4.9 K). IChapmaii et all l|2005l . = 36 ± 7 K), 
and lAmblard et all lj2010l . Td = 28 ± 8 K). The assump- 
tion /? = 2 follows the common practice in the literature, 
and is supported by many theoretical consid erations (e.g., 
iDraine fc Le3 1 19841 : iMathis fc WhiffenI Il989l l. The choice 
""i/m ~ 0-5 implies that the "efficient" mass range for star- 
formation covers about one and half orders of magnitude in 
halo mass. Though this number is somewhat arbitrary, it is 
unlikely to severely bias the final results, since, as demon- 
strated below, the power spectra are relatively insensitive to 
^2 

"L/m- 

We next vary one parameter at a time, and examine how 
the power spectrum changes compared to the fiducial case. 
The overall normalization Lq is fixed using the amplitude 
of the mean CIB measured by FIR AS (this assumption will 
be discussed and relaxed in ij7] below). More specifically, we 
choose Lo such that the ratio of the predicted and measured 
mean intensity, averaged over all frequency bands, is equal 
to unity. 



Nf 

1 ■^-^ IpiRAS.Vi 

iV7 ^ 1 

/ j — X model, 1^^ 



= 1. 



(26) 



Here Nf is the number of frequency bands in which the 
power spectra have been measured (for the most recent 
Planck measurement by All, Nf — 4), Imodei,^- is the 
intensity at frequency Vi predicted by our model while 
IpiB.AS,v is the same quantity measured by the FIRAS 
instrument. For the latte r ^ we use the values quoted in 
ICispert. Lagache fc Pugetl (|2000l l. 

The results are shown in Figure [1] for the highest (857 
GHz) and lowest (217 GHz) frequency bands of Planck. The 
power spectra in the fiducial model are plotted as solid thick 
curves, against which the other curves should be compared. 
There are two distinct ways in which the parameters can af- 
fect the power spectra: either by changing the shape of the 
spectrum (C^ oc 7^) or by changing the effective bias fac- 
tor (&). If the parameter enters primarily through the spec- 
trum, then it must have opposite effects at high and low 
frequencies, since the average intensity is kept fixed by the 
normalization (equation I26|l . On the other hand, if the pa- 
rameter enters primarily through the bias factor, then it has 
similar effects at different frequencies. This distinction helps 
us recognize how the parameters affect the power spectra, 
and also indicate their degeneracies - parameters which en- 
ter the same way are more likely to be degenerate. Figure [1] 
reveals that three of the parameters, Sz, Td and /3, affect 
the spectrum, while the other two parameters, rrie// and 
""i/mi enter through the bias factor. Consequently, degen- 
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Figure 1. CIB power spectra predicted by our models at 217 GHz (left panel) and 857 (right panel). The solid curve in both panels 
corresponds to our fiducial model, with parameter values of Sz = 3, me// = 10^'^-^ Mq, <j\i^ = 0.5, Td = 34 K, and /3 = 2 (see ^4Al 
for the definitions of these parameters). The other curves show how the power spectrum is modified when an individual parameter is 
changed by the amount indicated in the labels. 



eracies are likely strong among Sz, Td and /?, and between 
m,ff and a|,/^. 

Interestingly, increasing me// not only changes the 
overall clustering strength, but also increases the amplitude 
of the non-linear (1-halo term) contribution, relative to the 
linear (2-halo term) contribution to the power spectrum. 
This is seen clearly as the upturn in both panels of Figure [T] 
in the short-dashed [blue] curve at £ ^10'', and is similar to 
the effect of increasing the value of the parameter a in the 
previous models (see eguation llUp . In this sense, rrieff differs 
from Mmin , since the latter mostly affects the overall cluster- 
ing. This feature of our model arises from the non-flatness of 
the L — M relation, and is worth discussing further, because, 
as we will see below, the data requires significant power on 
small scales, which can be provided by the 1-halo term. 

For the intensity fluctuations, what matters for the rel- 
ative strength of the 1-halo and 2-halo terms is the total 
luminosity, rather than total number, of satellite galaxies. 
To be more concrete, let us suppose that the total lumi- 
nosity of satellite galaxies inside a halo of mass M can be 
parameterized by a power-law, similar to equation (|10|l . 



Ltot 



M 

Msat 



(27) 



It is then a_t, rather than a, that determines the small scale 
clustering. If the luminosities of galaxies were independent 
of the masses of their host (sub)halos, ul would be equal to 
Q. However, for a more realistic L — M relation, for which 
luminosity and mass are positively correlated over the rel- 
evant mass range, is larger than a. This is because in 
addition to the number of subhalos increasing with the halo 
mass, the average sub-halo mass increases, as well. Conse- 
quently, satellite galaxies are, on average, brighter in more 
massive halos, rendering ql larger than a. 



To illustrate this point explicitly, in Figure [2] we show 
the total number (solid curve) and the total luminosity 
(dashed and dotted curves) of satellite galaxies inside a 
halo of mass Mhaio- AH curves are shown normalized by 
their values at Mhaio = 10^^ Mq. As expected, the to- 
tal luminosity increases more rapidly with halo mass than 
the number of satellites. Further comparing the dashed 
(me// = 10^^-^ Mq) and the dotted (me// = 10^^ Mq) 
curves shows that the slope of the relation itself increases 
with rueff. This explains why the amplitude of the 1-halo 
term increases with me// compared to the 2-halo term. 

In their recent an alysis of the CIB power spectrum, 
lAmblard et al.l (|2011al ) have found that in order to fit the 
(large) observed power on small angular scales, the num- 
ber of satellite galaxies had to be increased compared to 
that expected from numerical simulations and optical sur- 
veys. In particular, the expectation is a <^1 and Msat = 
10 - 25M^» „ (e.g . . ICao et al .' 12004'; i Kravtsov et^ |2004| : 



IZheng et al. 20051 : iHansen et al. 2009), while the analysis 
of 'Ambl ard et al.f (f2011al ) requires either Msat — S.SiVfmin 
or, effectively, a > 1 (reducing Msat and increasing a have 
similar effects on the power spectrum since both increase 
the number of satellites and therefore raise the small scale 
power) . Our model naturally resolves this tension by distin- 
guishing a from ul, and, as we will demonstrate below, is 
able to fit the data without the presence of any additional 
low-mass satellite galaxies. 



5 CONSTRAINTS FROM PLANCK 

We are now ready to derive constraints from the Planck 
data. Given the flexibility of our extended model, we at- 
tempt to fit the CIB power spectra in the four Planck fre- 
quency bands simultaneously. As mentioned above, however. 
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Figure 3. CIB power spectra at 217 GHz, 353 GHz, 545 GHz, and 857 GHz measured by Planck (data points with error bars; taken 
from All) and predicted in our best— fit models (curves). The solid [red], long-dashed [green] and short-dashed [blue] curves show the 
total power spectra, including the shot noise, in case 0, 1 and 2, respectively (see Table 1 and ^for a summary of the models). For 
reference, the contributions from the shot noise (dotted curves), the 1-halo term (straight dash-dotted lines at £ > 10^) and the 2-halo 
term (dash-dotted curves) are shown separately in the fiducial "case 0" . 



Table 1. Best-fit values and marginalized l-a errors, as well 
as the reduced from fits to the Planck CIB power spectra. 
Within the parenthesis are the numbers of degrees of freedom. 
Five variants of our model are considered. Case 0: simultaneous 
fit to the 4 frequency bands of Planck, with three of the five model 
parameters allowed to vary, as explained in JS] Case 1: same as 
case 0, except s^ = enforced at z > 2. Case 2: same as case 
0, except the 857 GHz channels is ignored. Case 3: the graybody 
SED slope p is fixed at 2, and the dust temperature is instead 
a free parameter. Case 4-' the normalization Lq is allowed to vary 
(and found to be constrained to the range Lq = (0.73 ± 0.05)Li,, 
where Lt is the solution of eq. \26\ }. 





Sz 


logiorueff 


/3 


reduced 


Case 
Case 1 
Case 2 


4 64+°-" 
c c;7-l-0.99 
■^■"^'-0.96 
r, qQ+0.65 
'-'■^■^-0.65 


I2.65lg;l^ 
12.65lH^ 
12.78i°|0 


^■O»_0.08 

-"-■-"J-O.OT 
-, 7-1+0.15 


1.82 (33) 
2.78 (33) 
1.12 (24) 




Sz 


logiom^ff 




reduced x^ 


Case 3 
Case 4 


A oc-l-O.SO 
4-»0_o.79 

0-48_(, 44 


1 9 7O+0.23 
±^.(0_Q 23 

13.561^:^7 


44.8ll-?« 
43.49t^3:« 


2.43 (33) 
1.70 (32) 



there are strong degeneracies among our 5 model param- 
eters, which prevent us from simultaneously constraining 
each. We therefore first proceed by fixing two of the pa- 
rameters - the dust temperature Td, and the mass-range 



for significant luminosity, cr^y^ - at their fiducial values. In 
addition, we also apply a flat prior on s^, < Sz < 7. Here- 
after, this will be referred to as "case 0"; we will explore 
other possibilities below. 

To obtain the best-fit values and their confidence lev- 
els, we adapted the CosmoMC Mont e Carlo Markov Chain 
(MCMC) code l|Lewis fc Bridlell2002l ). and applied it to the 
data in the four Planck channels presented in All. The 
chains are checked visually and diagnosed using the statisti- 
cal tests provided by CosmoMC to ensure their convergence. 
The best-fit values and the marginalized errors are listed in 
the first row of Table [T] and the power spectrum in the best- 
fit model is shown explicitly at each of the four frequencies 
by the solid [red] curves in FigureO As this table and the fig- 
ure demonstrates, with our admittedly crude assumptions, 
and varying only three parameters, we are able to obtain 
good fits to the Planck data, with the exception of the high 
I data points at 545 GHz. where we underpredict the ob- 
served power. 

The best-fit Sz ~ 4.6 is qualitatively consistent with 
the rapid evolution of the sSFR indicated by recent measure- 
ments. As we will see below, the uncertainties induced by our 
model assumptions are likely much larger than the statistical 
errors. It is nevertheless interesting that the data, without 
any direct measurements of redshifts, requires a fairly rapid 
evolution of the sSFR. 

The halo mass scale for most efficient star formation 
is constrained to be around me// ~ 10^^'^^ Mq. This is 
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halo 



Figure 2. The solid curve shows the number of subhalos more 
massive than 10"'^"'^ M0 inside a halo of mass Mi^aiot while the 
dashed and dotted curves show the total luminosity of the these 
subhalos. All curves are normalized by their values at M^alo = 
10^2 M0, and m^ff is taken to be lO^^-^ Mq and lO" Mq for 
the dashed and dotted curves, respectively. The figure illustrates 
that the total luminosities of the satellites increase more rapidly 
than their number; this feature of our model helps to explain the 
large CIB power observed on small angular scales. 



consistent with the typical masses of the host halos of sub- 
millimetre galaxies, as predicted by semi-analytical models 
(e.g. Gonzalez et al. 2011). This mass scale is also in general 
agreement with that inferred from the clustering measure - 
ments of resolved bright sources (e.g., ICoorav et alll2010h . 
These studies all converge and indicate that star-formation 
is most vigorous in halos with masses of a few xlO^^ Mq. It 
is therefore interesting to note that a recent measurement of 
the small-scale CIB power by Herschel derive d a muc h lower 
mass-scale of 3 x lO" Mq. lAmblard et"all (|2011bl ) There 
are two possible reasons for this difference. First, the mea- 
surements themselves might be inconsistent. A recent com- 
parison between the two measurements indeed shows that 
the power spectra measured by Planck at 545 GHz and 857 
GHz are higher than those measured by Herschel. Second, as 
mentioned in the Introduction, the meaning of Mmin in the 
old model (adopted by the Herschel team) does not corre- 
spond to that of the rrieff defined here, invalidating a direct 
comparison between the two. 

The best-fit value for /3 (« 2.7) is larger than the value 
commonly assumed (/3 = 2). The common assumption is due 
to the fact that simple models for both insulating and con- 
ducting materials naturally give P = 2 at long wavelengths 



iDraine fc Leel 1 9841 : [Gordonll 19881 '). However, long; range dis- 
order in the dust grains can lead to /? > 2 (|Menv et al.l 
I2OO7I ). iMenv e t al. (200 J) also refers to observational and 
laboratory evidence for a wide range of /3 values, including 
those greater than 2. Note that others have also inferred a 
steep spectral index f or the background anisotropy spectrum 
at long wavelengths (|Hall et al.l 120101 : iDunklev et"alll2010l : 
IShirokofret"ai1l201ll ). 

Inferring (5 from observations is complicated by contri- 
butions from dust grains at varying temperatures. A mixture 
of cold and hot dust flattens the spectrum in the wavelength 
range between the two spectral peaks; if analyzed with a 



single-temperature grey-body model, the inferr ed value of P 
will t hus b e artificially low as emphasized by iReach et al.l 
(e.g.,|i99l). 

Just as a model that is missing a cold component can 
lead to an artificially low inference of /3, missing a hot com- 
ponent can lead to an artificially high inference of jS. What 
we infer as a high value of /3, may be due to a component, 
not included in our model, with a high effective temperature 
where Tcfi = Td/ (l + z). Perhaps we are missing a significant 
low-redshift component. 

More importantly, we emphasize that our best-fit /3 
value depends especially sensitively on model assumptions. 
In particular, as discussed above, /? is degenerate with s^; 
we will shortly see that /3 is reduced to ~ 2 if we make a 
different assumptions about the redshift-evolution. Another 
degeneracy of /3 is with Ta; a possible cause of a large /3 
is that we adopted the wrong dust temperature. The value 
Td = 34K is motivated by the temperature measurements in 
resolved sources, which, on average, are brighter, and reside 
at lower redshift, than the fainter sources responsible for the 
unresolved background. Finally, the measurement errors are 
considerable ( Jj5K). 

For the reasons in the preceding paragraph, we next 
consider the case where we fix /3 = 2, and allow the dust 
temperature to vary instead. The results are listed in Ta- 
ble [1] under "case 3". The values of s^ and logiorrisff are 
consistent with those from case 0, while the dust tempera- 
ture is found to be around 45 K, «11 K higher than assumed 
in case 0. This numerical result confirms a rough estimate 
based on Figure[T] a reduction of 0.68 in /3 could be compen- 
sated by an increase of ~ 10 K in T^. The best-fit dust tem- 
perature is then somewhat higher than the measurements 
quoted above (although within their 2cr errors). One could 
indeed speculate that at the higher redshifts probed by the 
CIB measurement, galaxies are more compact, and dust in 
their interstellar medium is hotter, since dust particles reside 
closer to the stars. 

In Figure |4l we show 2-D confidence contours of the 
parameters in case 0. These the correlations between the 3 
free parameters. In particular, rrieff is anti-correlated with 
Sz and P, while Sz and are positively correlated with each 
other. These correlations can be understood by examining 
the dependence of the power spectrum on individual param- 
eters, shown in Figure [T] above. 

Finally, owing to the current uncertainty on whether the 
sSFR continues to rise beyond redshift z > 2 (as discussed 
above) , we consider a variation on our fiducial model with a 
high-2 plateau. In this variant (referred to as "case 1"), we 
allow the L — M normalization to increase with redshift as 
before, i.e. as a power-law with a constant Sz feq|22[) at low 
redshift, but we set Sz = at z > 2. The results for this case 
are shown in the second row of Table [1] and by the dashed 
curves in Figure[3] The best-fit Sz and /3 both change consid- 
erably, with Sz increasing to Sz ~ 5.6, and j3 « 1.9 reduced 
close to its expected value of /3 ~ 2. This indicates that un- 
certainties in Sz and j3 due to model assumptions are large. 
Interestingly, the best-fit rrias remains almost unchanged. 
The minimum in case 1 is somewhat larger than in case 
(2.78 v.s. 1.82), indicating that the data favor case 0. 
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Figure 4. Constraints on tiie parameters Sz, mg.ff and /3 derived 
from Monte Carlo Markov Ctiains, by simultaneously fitting the 
CIB power spectrum in the four Planck frequency bands (in the 
fiducial model "case 0" defined in Q. The inner [red] and outer 
[yellow] contours show the 68% and 95% confidence levels, respec- 
tively. In each panel, the third parameter has been marginalized 
over. 



6 SHOT NOISE AND THE DUTY-CYCLE OF 
THE UNDERLYING INFRARED SOURCE 

As mentioned in i]3.2l above, our model can be used to ob- 
tain a lower limit on the shot noise. As we show here, this 
enables us to place a lower limit on the the duty-cycle of 
the underlying infrared sources, fduty, and to thus speculate 
on their nature (intermittent or quiescent). For simplicity, 
throughout the discussion below, we will assume a universal, 
constant duty cycle, independent of redshift and luminosity. 
Combining equations H19|) and (|20|) implies that 

0'"°* oc L^tNact, (28) 

where Lact is the luminosity of "active" galaxies, and Nact 
is the number of active galaxies. For a fixed abundance 
of galaxies, Nact is proportional to the duty-cycle, and if 
the total luminosity density (LactNact) is fixed, then Lact 
is inversely proportional to it. Thus, the amplitude of the 



shot noise power is inversely proportional to the duty-cycle, 

C"'"' oc 1/UutyE 

In Table [21 we list the shot noise estimates in our mod- 
els, in Case and Case 1 (the conclusions from the other 
cases are very similar) and compare these with results from 
the empirical model of Bll. The values of Scut needed in the 
calculation can be found in Table 3 of All. The Bll shot 
noise estimates have been carefully calibrated with source 
count measurements. The model employs a double expo- 
nential function to parametrize the luminosity function, and 
double power-laws to parameterize its redshift evolution. By 
adjusting 13 free parameters, the model successfully repro- 
duces the counts from the mid-infrared to the millimeter 
wavelengths, as well as the mid-infrared luminosity func- 
tions. The duty cycle is taken to be either fduty = 1 (mim- 
icking the value for quiescent star-formation in long-lived 
galaxies) or fduty = 0.01 (mimicking short-lived star-bursts, 
as expected when star-formation is triggered by major merg- 
ers). As the table shows, the shot noise levels for fduty ~ 1 
are of the same order as the results of Bll (although our 
shot noise levels for Case 0(1) are somewhat higher(lower) 
than those of Bll). On the other hand, with fduty ~ 0.01, 
the shot noise in our models severely exceed those of Bll 
(by factors of ^20). Given this vast difference, we con- 
clude that the average duty cycle of the underlying sources 
is on the order of unity. This suggests that the infrared back- 
ground is mainly contributed by normal quiescent galaxies, 
rather than starbursts. This con clusion is consistent with re- 
cent numerical simulations (e.g.. lHopkins et al.ll20ld ). which 
find that ~ 85% of the infrared background is contributed 
by normal galaxies. 

It is interesting to note that between the two cases, our 
Case 1 has a lower shot noise, due to its flatter redshift evolu- 
tion. With the total intensity fixed, the shot noise decreases 
with increasing source number density. In Case 1, a higher 
fraction of the CIB originates from low redshift, where the 
galaxy number density is higher. 

It might appear worrisome that our lower limits on the 
shot noise, in Case 0, are slightly higher than those of Bll. 
This, however, could arise because our assumed (t^^^ is too 
small. Increasing cr| would reduce the average luminosity, 
and consequently reduce the shot noise levels. If this were the 
case, then me// in Table[T]would have been underestimated. 
On the other hand, comparisons with the BLAST measure- 
ment show that the shot noise estimates by the model of 
Bll could be lower than the true values by a considerable 
amount (60% at 545 GHz, 20% at 857 GHz; see All). 

In principle, our IR source population model could be 
compared directly to measurements of the counts of resolved 
bright sources, as well. The distribution of source counts as 
a function of flux, especially at the bright end, also depends 
sensitively on the assumed the duty cycle. In particular, if 
the total CIB intensity is kept fixed, then the number of 
bright sources increases, while the number of faint sources 



° This neglects the fact that the brightest sources have been 
masked. For Planck, the fiux cut is large (a few hundred mjy), 
and we find that the inverse relation between C'"''"* and fduty 
holds accurately as long as fduty > 0.1. 
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Table 2. Shot noise levels obtained by Bll, compared with the 
lower limits computed in our models, for Cases and 1, and for 
two different values of the duty-cycle fduty 





fduty 


217 GHz 


353 GHz 


545 GHz 


857 GHz 


Bfl 




12.2 ±2.9 


138 ± 22 


1150 ±92 


5923 ± 367 


Case 


1 

0.01 


14.79 
516.93 


308.39 
5790.21 


2305.94 
34813.30 


7874.02 
85269.09 


Case 1 


1 


5.71 


114.43 


1057.59 


5434.20 


0.01 


570.56 


7884.97 


72353.10 


89765.42 



decreases with decreasing duty-cyclel3 Because we do not 
include any scatter in the L — M relation, which would in- 
crease the number of bright sources, the source counts pre- 
dicted by our model are lower limits at the bright end. By 
comparing the model predictions with measurements, we are 
therefore able to derive a lower limit on the duty-cycle. We 
show the results of such a comparison at 857 GHz in Fig- 
ure (5] The data points are taken from the Herschel Mult i- 
tiered Extra-galactic (HerMES) survey (|01iver et al .l2010bh . 
and the long-dashed, short-dashed, and dotted curves are 
the lower limits on the counts, computed in our model for 
fduty = 1, 0.1 and 0.01, respectively (with the other param- 
eters fixed at their best-fit values in Case 0). This figure 
clearly shows that only fduty = 1 is consistent with the 
observed number counts; the other two values significantly 
overpredict the number of the bright sources (by more than 
an order of magnitude) . We therefore arrive at the same con- 
clusion as that drawn from the analysis of the shot noise, 
namely that the CIB must be produced primarily by long- 
lived sources with fduty ~ 1- 



7 DISCUSSION 

In this section, we discuss a few remaining issues and caveats 
related to our main results. 

(1) Small scale clustering at 545 GHz 

Although we are able to fit all four Planck channels with 
an overall reduced of 1.82, formally this is still a bad fit, 
and our best-fit model noticeably underpredicts the small- 
scale power at 545 GHz (by a factor of ~two). There are at 
least two possible ways to alleviate this problem. 

First, the shot noise, which is almost parallel to the 
1-halo term over the scales considered in this study, is un- 
certain. The model could match the data better if the shot 
noise were raised by a suitable amount. Indeed, as we men- 
tioned previously, the shot noise levels measured by BLAST 
l|Viero et al.ll2009l ') are about 60% and 20% higher than the 
model predictions by Bll at 545 GHz and 857 GHz, respec- 
tively (All). If we keep the best-fit model parameters in 
Case as listed in the first row of Table [T] but increase the 
shot noise at 545 GHz and 857 GHz by these factors, the 



^ This can be easily verified by visualizing shifting a logAf-logS 
diagram "downward" and then toward higher luminosities. 
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Figure 5. Resolved bright source counts at 857 GHz from the 
HerMES survey (data points) are compared with lower limits de- 
rived from our models, assuming fduty = 1 (long dashed [green] 
curve), 0.1 (short-dashed [blue] curve) and 0.01 (dotted [red] 
curve) . The other model parameters are set to their best-fit values 
in Case 0. 



tension between the model and data at the small scales is 
significantly reduced. 

Second, one can certainly attempt to modify and/or 
fine-tune our best-fit models, to increase the small-scale 
power. As we mentioned in SQ] this could be achieved by 
increasing rrieff. However, we have found no simple way 
to increase the small scale clustering at 545 GHz, without 
affecting other frequencies, especially at 857 GHz. In partic- 
ular, this is because there is considerable overlap between 
different frequencies, caused by the extended range of red- 
shifts over which sources contribute to the CIB. To illustrate 
this inter-dependence of frequencies, we perform a fit to the 
three lowest frequency bands, ignoring the 857GHz channel 
(denoted as "Case 2" in the third row of Table [TJ and shown 
with short dashed curves in Figure The best fit to these 
three channels yields a oi 27 for 24 degrees of freedom. 
The power at 857 GHz in this case is, however, overpredicted 
by a factor of Ritwo. We conclude that non-trivial modifica- 
tions to the SED, or other ingredients of our model would 
be needed in order to improve this situation. 

(2) The normalization of the L — M relation 

Thus far, the normalization of the L — M has been fixed 
by equation (|26p . This makes the MCMC chains converge 
more easily, and yields tighter constraints on other param- 
eters. However, fixing Lo with infinite precision could raise 
two possible concerns. First, while equation (|26|l guarantees 
that our model is on average consistent with the existing FI- 
RAS measurement of the mean CIB, individual frequencies 
might still deviate, with a few frequencies far below the mea- 
sured values, while the others far above. Second, since the 
uncertainties of the F77?^>S measurement are non-negligible, 
could our conclusions be altered if the normalization is al- 
lowed to vary? 

To address the first concern, in Figure |6] we show the 
spectrum of the mean CIB, as predicted in our model 
(for Case only, as the spectra in the other cases are 
very similar). This figure also shows the FIR AS measure- 
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Figure 6. The spectrum of the mean CIB in our fiducial Case 
0, shown at the four frequency bands of Planck (sohd squares). 
The spectra in our other models is similar. Our mean background 
is consistent with the background spectrum measured by FIR AS 
(solid curve) , lying w ell within the ±lcr a llowed region (indicated 
by the dotted curves: [Cispert. Lagache fc Pugelii2000l ). 



ment, along with the ilcr region (as derived from the 
errors on the parameters in t he FIR AS fitting formula; 
iGispert. Lagache fc Puge3l200d ). Clearly, the intensities at 
all frequencies are consistent with the observations. 

To address the second concern, we relaxed the normal- 
ization Lo in the Monte-Carlo fitting. Other settings were 
kept the same as in Case 3, except that we imposed the prior 
nieff ^ 10^* M0, in addition to the prior on Sz- This was 
necessary because otherwise the fits allowed combinations 
of unphysically high m^f f with suitably low normalizations. 
With the addition of Lq, we then have a total of four free 
parameters. Only the Planck data was used in the fitting, 
although in principle, the FIRAS data could be used explic- 
itly for the normalization, as well. The marginalized errors 
on the parameters are shown in Table [T] marked "Case 4". 
The best-fit value for the normalization Lq is found to be 
(0.73 ± 0.05)1/*, where L* is the normalization satisfying 
equation (|26|) . Compared to other cases, rUeff ~ lO'^"^'® 
is considerably higher and well beyond the halo mass of a 
typical galaxy. The large me// raises the small scale clus- 
tering, and helps alleviate the small scale tension we dis- 
cussed previously. It also increases the large-scale cluster- 
ing (by a smaller factor), which is then compensated by 
the lower normalization Lq. The negative correlation be- 
tween me/ / and Lq is also clear from the confidence contours 
shown in Figure [T] Interestingly, the best-fit normalization, 
~ 0.73 Ifiras, is within the allowed range of the FIRAS 
measurement, despite the fact that this was not imposed on 
the fits. 

(3) Future improvements of the model 

Although successfully describing the Planck data, and 
an improvement over previous version, our model is still ad- 
mittedly very simplified. In the following, we give a few ex- 
amples how the model could be improved and extended to 
include more physical ingredients. 

1) As mentioned above, we neglect any scatter in the 
L — M relation. This scatter could be modeled, and would 



Figure 7. Confidence contours (68% and 95%) on logxom^ff 
and Lq for Case 4, in which the normalization Lq was allowed 
to vary. The contours have been marginalized over the other two 
free parameters, Sz and T^. 



impact the estimates of the shot noise and the duty-cycle, in 
particular. The scatter itself might be mass- and redshift- 
dependent. 

2) We have a single galaxy population in our model. In 
reality, multiple populations, such as normal blue galaxies, 
star-bursts and obscured AGNs may all make appreciable 
contributions to the CIB. These populations might have dif- 
ferent SEDs, duty cycles, L — M relations and redshift evo- 
lutions. The fraction of starbursts and AGNs could be taken 
from theoretical merger-tree models, calibrated by existing 
data on the IR galaxy populations. In addition, the vari- 
ous properties of each population may depend not only on 
redshift, but also on environment. 

3) In the current model, very simple forms have been 
assumed for the SED and L — M relations, and the SED was 
not allowed to evolve with redshift. These are features that 
could be improved in future work. For example, an (asym- 
metric) double power-law with a suitable set of parameters 
is probably a better representative of the true L — M rela- 
tion, though it requires more parameters. 

4) Central and satellite galaxies could have been treated 
with difi'erent L — M relations. 

With such improvements, the model could be developed 
so as to properly compute the shot noise, and simultaneously 
fit number counts of resolved sources and the power spectra 
of background fluctuations. 



8 CONCLUSIONS 

In this study, we developed a halo-model based formalism 
to compute the power spectrum of the cosmic infrared back- 
ground (CIB). Although previous, similar, models provided 
excellent fits to individual frequency bands, they were based 
on an unrealistic assumption that the luminosities of galax- 
ies are independent of the host halo masses. We relaxed this 
assumption by incorporating the subhalo mass function, to- 
gether with a more realistic L — M relation, into our model. 
With these improvements, we were able to naturally resolve 
the tension between the high power in the CIB on small an- 
gular scales (requiring a > 1, where the number of satellite 
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galaxies A*'sat in a halo scales as oc M") and optical data and 
numerical simulations results (which favor a <,!)■ 

With our improved model, we were also able to simulta- 
neously fit all four frequency bands of the Planck measure- 
ment of the CIB anisotropy spectrum, by varying only three 
parameters: Sz (the describing the redshift evolution of the 
L — M relation), rUeff (the halo mass-scale on which star- 
formation is most efficient), and P (a parameter describing 
the graybody SED of the unresolved IR sources) . We found 
that star formation is most efficient in halos with relatively 
high mcisses of rUeff =severalxlO^^ M0, in general agree- 
ment with semi-analytical models of galaxy formation. We 
also found that the Planck data favor an increase in the 
L — M normalization with redshift. Finally, by comparing 
the shot noise and bright source counts predicted by our 
models with an observationally calibrated empirical model 
and a direct measurement, respectively, we conclude that 
the duty-cycle of the source is on the order of unity. This 
implies that the CIB is dominated by long-lived, quiescent 
star-forming galaxies, rather than short-lived star-bursts. 
This conclusion should be robust, since our estimates of shot 
noise and bright source counts are already conservative, ow- 
ing to our neglect of any scatter in the L — M relation. 

The model we presented here, with the potential to com- 
pute the shot noise and number counts, in addition to the 
correlated background fluctuations, provides a theoretical 
framework for a future joint analysis of the unresolved CIB 
background and resolved source counts. Such a coherent 
analysis will no doubt afford new insight into the relation 
between dust-enshrouded star-formation and dark matter 
halo properties. 
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